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We analyze the parametric space of the constrained minimal supersymmetric standard model with 
/i > 0 supplemented by a generalized asymptotic Yukawa coupling quasi-unification condition which 
yields acceptable masses for the fermions of the third family. We impose constraints from the cold 
dark matter abundance in the universe and its direct detection experiments, the B-physics, as well 
as the masses of the sparticles and the lightest neutral CP-even Higgs boson. Fixing the mass of 
the latter to its central value from the LHC and taking 40 < tan /3 < 50, we find a relatively wide 
allowed parameter space with —11 < AolMxji < 15 and mass of the lightest sparticle in the range 
(0.09 — 1.1) TeV. This sparticle is possibly detectable by the present cold dark matter direct search 
experiments. The required fine-tuning for the electroweak symmetry breaking is much milder than 
the one needed in the neutralino-stau coannihilation region of the same model. 
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I. INTRODUCTION 

The well-known constrained minimal supersymmetric 
standard model (CMSSM) [J-Q is a highly predictive 
version of the minimal supersymmetric standard model 
(MSSM). Its basic characteristic is that it employs uni¬ 
versal boundary conditions for the soft supersymmetry 
(SUSY) breaking terms. Namely, the free parameters of 
the CMSSM are the following 

sign^, tan/3, M 1 / 2 , mo, and Aq, (1) 

where sign/i is the sign of /i, the mass parameter mix¬ 
ing the electroweak Higgs superfields H 2 and Hi of the 
MSSM which couple to the up- and down-type quarks 
respectively, tan /3 is the ratio of the vacuum expecta¬ 
tion values of H 2 and Hi, while the remaining symbols 
denote the common gaugino mass, the common scalar 
mass, and the common trilinear scalar coupling constant, 
respectively, defined at the grand unified theory (GUT) 
scale Mgut, which is determined by the unification of 
the gauge coupling constants. 

However, the recently announced experimental data 
on the mass of the standard model (SM)-like Higgs bo¬ 
son and the branching ratio HR (Hg —>■ /i+p, ) of 

the process B„ [lll - H^ in conjunction with the 

upper bound [l4| from cold dark matter (CDM) consid¬ 
erations on the relic abundance of the lightest 

SUSY particle (LSP) - which is the lightest neutralino 
X - put under considerable stress (la - H^ the parameter 
space of the GMSSM. It is by now clear, however, that 
there still exist viable slices of the following three regions 
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of the parameter space of the CMSSM which give accept- 
able values of Hlsp^^: 

(i) The neutralino-stau (x — ti) coannihilation region 
with Mi /2 ^ mo, where y is an almost pure bino 
and Hnsp/i^ is reduced to an acceptable level due to 
the proximity between the mass mnsp of the LSP 
and the mass mf^ of the lightest stau fi, which is 
the next-to-LSP (NLSP). 

(ii) The A-funnel region [l^ appearing at large (> 40) 
tan/3 values, where the mass m .4 of the GP-odd 
Higgs boson A satisfies the relation niA — 2?71lsp 
and so the y — y pair annihilation procedure is en¬ 
hanced by an H-pole exchange in the s-channel, 
thereby reducing Hlsp^^- 

(iii) The hyperbolic branch/focus point (HB/FP) region 

|19i[27l | located at large values of mo Mi /2 and 

small /x’s, which ensures a sizable higgsino frac¬ 
tion of y. The Hlsp^^ remains under control for 
mpsp ^ 1 TeV thanks to the rapid y ~ X annihila¬ 
tion and the neutralino-chargino (y/y 2 ~yi') coan¬ 
nihilation [ 2 § - l^ (y 2 is the next-to-lightest neu¬ 
tralino and the lightest chargino). 

It would be interesting to investigate whether this 
available parameter space survives in even more restric¬ 
tive versions of the CMSSM, which can emerge by em¬ 
bedding it in concrete SUSY GUT models. Here we focus 
on the Pati-Salam (PS) GUT model based on the gauge 
poup Gps = SU{4)c X S't/(2)L X SU{2)r P, Ipl. It 
is interesting to note that this model can arise [33, |34| 
from the standard weakly coupled heterotic string and 
the D-brane setups P, P] and its phenomenology has 
been extensively studied - see e.g. Ref. [p. It can also 
naturall y a rise in the framework of non-commutative ge¬ 
ometry [381 . 
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The SUSY PS model in its simplest realization leads 
to ‘asymptotic’ Yukawa unification (YU) [13, 1^ . 
i.e. the exact unification, at the GUT scale Mgut, of 
the Yukawa coupling constants ht, ht, and hr of the top 
quark t, the bottom quark b, and the tau lepton r, re¬ 
spectively. Although the CMSSM with YU is an elegant 
model, it yields unacceptable values of the 6 -quark mass 
mb for both signs of ^ given the experimental values of 
the t and r masses - which, combined with YU, naturally 
restrict tan (3 to large values. This is due to the presence 
of sizable SUSY corrections [42l - l45j to the tree-level b- 
quark mass (about 20 %), which drive it well above [a 
little below] its 95% confidence level (c.l.) experimental 
range for ij, > 0 [/r < 0 ] - we use the standard convention 
of Ref. for the sign of fx. 

In order to circumvent this difficulty, we consider 
here the extended PS SUSY GUT model introduced in 
Ref. [13 and reviewed in Refs. [H, ji^, which yields a 
moderate deviation from exact YU and, thus, allows ac¬ 
ceptable values of the 6 -quark mass for both signs of 
within the GMSSM. In particular, the Higgs sector of 
the simplest PS model [3l|, jsll is extended so that H 2 
and Hi are not exclusively contained in a S'?7(4)c sin¬ 
glet, S'?7(2)l X S'C/(2)r bidoublet superfield, but receive 
subdominant contributions from another bidoublet too 
which belongs to the adjoint representation of S'C/(4)c. 
As a result, YU is naturally violated and replaced by a 
set of asymptotic Yukawa quasi-unification (YQU) con¬ 
ditions: 


/it(MGUT) : hbiMaur) ■ hriMcvr) = 


1 - pa2l'Ji 


I — pail\fZ 


1 - 1 - \fZpai 

v'l -f |a 2 P 


v'l + laiP 


•\/l + |cnP 


These conditions depend on two complex parameters ai, 
a 2 and one real and positive parameter p. The parame¬ 
ters ai and a 2 describe the mixing of the components of 
the SU{4:)c singlet and 15-plet Higgs bidoublets, while p 
is the ratio of their respective Yukawa coupling constants 
to the fermions of the third fam ily j^ . lb (i|. 

In our original papers [43, Isil - I^ , we have considered 
monoparametric versions of the YQU conditions depend¬ 
ing on just one parameter, which was considered for sim¬ 
plicity real and replaced tan fi in the investigation of the 
parametric space of the model. These forms of the YQU 
conditions arise by taking ai = —02 for /r > 0 or ai = a 2 
for p < Q. Indeed, these choices turn out to be suitable 
for generating an adequate violation of YU and ensur¬ 
ing, at the same time, successful radiative electroweak 
symmetry breaking (EWSB) and the presence of a neu- 
tralino LSP in a large fraction of the parametric space. 
However, the emergent versions of GMSSM are by now 
experimentally excluded. For p < 0, this is due to the in¬ 
compatibility between the bound on the branching ratio 
HR (6 —^ sy) of the process 6 —57 and the CDM con¬ 
straints [^. For /r > 0, on the other hand, it is due to the 
fact that the parameter space where the GDM constraint 
on is satisfied thanks to x — ti coannihilations 


turns out to be non-overlapping with the space allowed 
by the data [ni[i3 on BR(Hs ) (s^ j. The main 

reason for the latter negative result is that tan fi remains 
large and, hence, it enhances the SUSY contribution to 
BR{Bs -)> p~''p~)- 

In order to overcome this hurdle, we adopted in 
Refs. [4l,[53 the more general version of Eq. (I2|) without 
imposing any relation between ai and a 2 - As a conse¬ 
quence, we could accommodate more general values of 
the ratios hm/hn with m,n = t,b,T - still of order unity 
for natural values of the model parameters - and lower 
tan P’s. The allowed parameter space of the model found 
in Ref. [^ is then mainly determined by the interplay be¬ 
tween the constraints from the BR(i?s —>■ p '^ p ~), CDM, 
and the results of LHC on the Higgs boson h mass mh- 
Actually, our acceptable solutions are due to a synergy 
between x — ti coannihilation and the A-funnel mecha¬ 
nism - see points (i) and (ii) above - given that mn — mA 
{H is the heavier CP-even neutral Higgs boson with mass 
mu). However, tolsp comes out to be large (^ 1 TeV), 
which makes the direct detectability of the LSP very diffi¬ 
cult and the sparticle spectrum very heavy. Furthermore, 
the emergent values of M 1 / 2 , mo, and p lie in the multi- 
TeV range, which puts under some stress the naturalness 
of the radiative EWSB. 

Therefore, the tantalizing question, which we address 
here, is whether the conditions in Eq. @ can become 
compatible with the solutions obtained in the HB/FP re¬ 
gion of the CMSSM - see point (iii) above. We show that 
this is not only possible but also quite promising since, 
apart from allowing a much wider parametric space, it 
gives smaller values of p and wlsP) despite the fact that 
mo and M 1/2 are quite large and the sfermions, sleptons 
and heavier Higgs bosons are quite heavy. Moreover, the 
LSP can be probed by the CDM detection experiments, 
which play a crucial role in limiting the parameter space 
of the model. Finally, the deviation needed from exact 
YU can be naturally attributed to the YQU conditions 
in Eq. while the fine-tuning required for the radia¬ 
tive EWSB turns out to be milder than the one needed 
in the X — ti coannihilation region for musp’s yielding 
Hlsp^^’s close to the cosmological upper bound. 

We first (Sec.ini review the phenomenological and cos¬ 
mological requirements which we consider in our investi¬ 
gation. Next (Sec. IHII) . we exhibit the salient features of 
the HB/FP region of the CMSSM and find, in Sec. lIVl the 
resulting restrictions on the parameters of the CMSSM. 
We, finally, check, in Sec.|Vl the consistency with Eq. ([2]) 
and discuss the naturalness of the model in Sec. ED We 
summarize our conclusions in Sec. EID 


II. PHENOMENOLOGICAL AND 
COSMOLOGICAL CONSTRAINTS 

In our investigation, we closely follow the renormaliza¬ 
tion group and radiative EWSB analysis of Refs. [4^.[50j. 
Most notably, we use an optimal, common supersymmet- 
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ric threshold 


2. The mass of the lightest Higgs boson 


Msvsy ^ 


( 3 ) 


(with ti 2 being the stop quark mass eigenstates) which 
reduces drastically the one-loop corrections to the Higgs 
boson masses, making their calculation from the tree- 
level effective potential quite accurate and stable [ 3 ] ■ The 
SUSY spectrum is also evaluated at Msusy by using the 
publicly available calculator SOFTSUSY M|. The output 
is put into micrOMEGAs a publicly available code, 
which calculates a number of phenomenological - see 
Sec. Ill Al - and cosmological - see Sec. IIIBI - observ¬ 
ables assisting us to restrict the parametric space of our 
model. 


A. Phenomenological Constraints 

The phenomenological considerations which we take 
into account in our study are the following: 


1. SM fermion masses 

In order to determine with good precision the running 
of the Yukawa coupling constants, we have to take into 
account the sizable SUSY corrections [s^, to 

the tree-level 5-quark and r-lepton masses. We incorpo¬ 
rate them in our code using the formulas of Ref. in 
accordance with the recent reanalysis of Ref. [1^. The 
result is to be compared with the experimental value of 
mb{Mz)- This is derived by using as an input parameter 
the MS value, which, at 1 — cr, is 

=4.18 ±0.03 GeV. (4) 

This range is evolved up to Mz using the central value 
as{Mz) = 0.1185 [1^ of the strong fine-structure con¬ 
stant at Mz and then converted to the DR scheme in 
accordance with the analysis of Ref. [13 . We obtain 

2.8 <mb{Mz)/GeY <2.86 (5) 

at 68% c.l. with the central value being mb{Mz) = 
2.83 GeV. Less important but not negligible (almost 4%) 
are the SUSY corrections [131 to the r-lepton mass, which 
have the effect [43,[ll| of slightly reducing tan f3. We take 
as input the central value |60| 

mr{Mz) = 1.748 GeV (6) 

of the DR r-lepton mass at Mz. Finally, as regards the 
t-quark mass, we take the latest 1 — ct result [6l| on its 
pole mass [Mt = 173.34 ± 0.76 GeV) and construct the 
68 % c.l. range for its running mass 

164.1 < mtimt)/GeV < 165.56 (7) 


The experiments ATLAS and GMS in the LHG discov¬ 
ered simultaneously a boson that looks very much like the 
expected SM Higgs boson. Its measured mass turns out 
to be [3, [l3| 

_ f 125.36 ± 0.37 (stat.) ± 0.18 (sys.) GeV (ATLAS), 
“ |125.03;°J? (stat.)t°;i3 (sys.) GeV (GMS). 

( 8 ) 

While there is no combined analysis of the two experi¬ 
ments yet, we can estimate the allowed 95% c.l. range 
of this mass including a theoretical uncertainty of about 
±1.5 GeV. This gives 

122 <TO/,/GeV< 128.5. (9) 

In most of the plots in Sec. El we set 

mb = 125.5 GeV, (10) 

which lies well inside the experimental range from both 
ATLAS and GMS in Eq. dH). 


3. Sparticles searches 

The mass limits which are relevant for our investigation 
are the following: 

a. Mass of the chargino The combined results of ex¬ 
periments that took place at LEP showed that, re¬ 
gardless of the model, the lower bound on the mass of 
the charginos is 

m^± > 103.5 GeV. (11) 

b. Mass of the gluino The present model dependent 
95% c.l. lower bound on the mass of the gluino is [ill 

m-g > 1.3 TeV. (12) 

It is well-known that such a heavy gluino spoils the suc¬ 
cess of several models which consider the HB/EP region 
[lUnl- However, as we will show, this does not happen 
in our case. 


f. B-physics constraints 


We also consider the following constraints originating 
from H-meson physics: 

• The branching ratio HR (Hg — of th® Pro¬ 
cess Bs -T ht~ - [zlizi is to be consistent with 

the 95% c.l. experimental upper bound [13, [T3| 

BR(Hs ^/i+/r-) < 4.2 X 10”®. (13) 


Although the most recent experimental result is 


with the central value being mt{mt) = 164.83 GeV. 


1.1 < BR(Hs ^ ^+/i-)/10-® < 6.4, (14) 
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we use the more stringent upper bound in Eq. 0, 
since it is a combined result and, thus, much more 
reliable. 

• The branching ratio BR (6 —>■ sod of & —>■ S7 |ljH73 
must lie in the 95% c.l. range [78 - [80l| 

2.79 X 10"^ < BR (6 ^ sq) < 4.07 x IQ-^. (15) 

• The ratio R (i?„ —)> tv) of the CMSSM to the SM 
branching ratio of the process Bu —>■ t v 
should be confined in the 95% c.l. range 

0.52 < R(B„ ^ Ti/) < 2.04. (16) 

5. Muon anomalous magnetic moment 

There is a discrepancy 5a^ between the measured value 
of the muon anomalous magnetic moment from its 
SM prediction. The latter, though, is not yet completely 
stabilized mainly because of the ambiguities in the cal¬ 
culation of the hadronic vacuum-polarization contribu¬ 
tion. According to the evaluation of this contribution in 
Ref. [13 , there is still a discrepancy between the findings 
based on the e^e“-annihilation data and the ones based 
on the r-decay data. However, in Ref. [s^, it is claimed 
that this discrepancy can be alleviated. Taking into ac¬ 
count the recent and more reliable calculation based on 
the e+e“ data [s^, the complete tenth-order QED con¬ 
tribution [s^, and the experimental measurements [s^ 
of a^, we end up with a 2.9 — cr discrepancy 

5a^ = (24.9 ± 8.7) x 10"^°, (17a) 

resulting to the following 95% c.l. range: 

7.5 X 10"^° < 5a^ < 42.3 x 10"^°. (17b) 

This 5a^ can be attributed to SUSY contributions cal¬ 
culated by using the formulas of Ref. [s^ . The resulting 
5a^ has the sign of /x and its absolute value decreases as 
wlsp increases. Therefore, Eq. (I17bl) hints that the sign 
of fi has to be positive. Moreover, a lower [upper] bound 
on mnsp can be derived from the upper [lower] bound 
in Eq. (Il7b|) . As it turns out, only the upper bound on 
wlsp is relevant here. Taking into account the afore¬ 
mentioned computational instabilities and the fact that 
a discrepancy at the level of about 3 — a cannot firmly 
establish a real deviation from the SM value, we restrict 
ourselves to just mentioning at which level Eq. (I17al) is 
satisfied in the parameter space allowed by all the other 
constraints - cf. Refs. |15l-[l^. 

B. Cosmological Constraints 

Our cosmological considerations include information 
from the CDM abundance and direct detection experi¬ 
ments quoted below. 




FIG. 1: Dominant x—X S'^inihilation and x/x 2 — Xi' coannihilation 

reactions in the HB/FP region. 

1. Cold dark matter abundance 

In accordance with the recently reported results [13 ] 
from the Planck satellite, the 95% c.l. range for the CDM 
abundance in the universe is 

= 0.1199 ± 0.0054. (18) 

In the context of the CMSSM, the lightest neutralino x 
can be the LSP and, thus, naturally arises as a CDM 
candidate as long as its relic abundance does not exceed 
the upper bound in Eq. 0 , i.e. 

^■Lsph"^ < 0.125. (19) 

This is a quite strong restriction on the parameter space, 
since DLSph^ increases, in general, with mLgp and so 
an upper bound on mpsp can be derived from Eq. 0 . 
Note that no lower bound on Dnsph^ is imposed in our 
analysis, since the CDM may receive contributions from 
other particles too 

Focusing on the HB/FP region, we find that the main 
processes causing the reduction of Onsph^ are the x~ X 
annihilations (for large \Aq/Mi/ 2 \^s and low mpsp’s) and 
the x/X 2 — xt coannihilations (for low \Ao/Mi/ 2 \'s and 
large mpsp’s). Specifically, the dominant reactions are 
(see Fig.[T|) - the notation used for the various (s)particles 
is the same as in Table UTl - 

• XX W~^W~ mediated by a /i and H exchange in 
the s-channel and 

• x/X 2 Xi" and ud mediated by a 61^2 or dL,R 

exchange in the t-channel, respectively. 

Here, as in Sec. lIVDi we consider the squarks of 
the two lightest generations as degenerate. The strength 
of the coannihilation processes is controlled by the rela¬ 
tive mass splittings or Ax 2 between the x^ or X 2 , 
respectively, and the LSP, which are defined by 

Aca = (wcA - wlsp)/wlsp, CA = x^,x2- (20) 

In particular, we find that the resulting Dnsph^ decreases 
with Ax± or A^s. The x/X 2 “ Xi" coannihilation pro¬ 
cesses are dominant in the regions where the bound on 
f^LSP^^ in Eq. (fTOl) is approached contributing almost 
40% to the total effective cross section. On the other 
hand, in the region where m-x* near its lower limit in 
Eq. (HID , the x~X annihilation processes contribute more 
than 50% to the reduction of Dpsph^ - cf. Ref. [s^. 
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2. Spin-independent direct detection of CDM 


Direct detection of dark matter through its elastic scat¬ 
tering off atomic nuclei inside an underground detector 
would be an undeniable evidence of its existence. Data 
coming from direct detection experiments can provide 
strict bounds on the values of the free parameters in 
SUSY models with a sizable higgsino component, as in 
our case. Here, we focus on the large underground Xenon 
(LUX) experiment 93|, whose present data are a little 


more restrictive than those from the XENONIOO exper¬ 
iment [ 9 ^. Both belongs to the experiments probing for 
the spin-independent (SI) cross section, since they con¬ 
sider only the scalar part of the cross section. Therefore, 
the quantity which has to be considered for comparing 
the experimental results and theoretical predictions is the 
SI neutralino-proton (y — p) elastic scattering cross sec¬ 
tion cr|p. This quantity is calculated by employing the 
relevant routine of the micrOMEGAs package [9^ based 
on the full one-loop treatment of Ref. [9^. The scalar 
form factors (with q = u,d, s) for light quarks in the 
proton, which are needed for the calculation of are 
estimated following the revised approach of Ref. [ 94 I and 
taking into account the recent lattice simulations [^ . 
which suggest the following 68% c.l. ranges for the pion- 
nucleon and the strangeness-nucleon {os) sigma 

terms: 


a^N = 45 ± 6 MeV and a, = 21 ±6 MeV. (21) 

Using also as inputs the light quark mass ratios 
niu/nid = 0.48 ± 0.1 and ms/md = 19.5 ± 2.5, we find 


the following I — cr ranges for the 


fP _ (Y niQ+0.0057 

JTu — ^•'^'‘-^-0.00445 

(22a) 

fP _ n Q25+0.0016 

(22b) 

= 0.022 ± 0.0064. 

(22c) 


Note that these ranges are much narrower than the ones 
used in Ref. [s^. Also turns out to be considerably 
smaller than its older value - cf. Ref. [ 2 ^ - reducing 
thereby the extracted cr|p. 

In the HB/FP region, where the LSP has a signifi¬ 
cant higgsino component, cr|p is dominated by the t- 
channel Higgs-boson-exchange diagram contributing to 
the neutralino-quark elastic scattering process - for the 
relevant tree-level interaction terms see e.g. the appendix 
of Ref. [ 4 ^. Especially for large tan/3’s, which is the 
case here, the couplings of the heavier Higgs boson H to 
down-type quarks are proportional to tan/3 and so are 
the dominant ones. More explicitly, cr|p behaves as 

cr~p oc tan^/3|Ai4p|Ni,3p/m^, (23) 

where IVi.i, IVi, 2 , and IVi,( 3 , 4 ) are the elements of the 
matrix N which diagonalizes the neutralino mass matrix 
and express the bino, wino, and higgsino component of 


X, respectively. As a consequence, <7-^ can be rather en¬ 
hanced compared to its value in the x — ti coannihilation 
region, where x is essentially a pure bino. 

The data from the LUX experiment, which we take 
from Ref. [9^, are directly applicable in the case where 
the CDM consists of just one kind of weakly interacting 
massive particle, which in our case would be the neu¬ 
tralino X- However, if the x’s constitute only a part of 
the CDM in the universe, some extra ca re is requ ired. 
The number of the scattering events A is I0Cll - ll02j| pro¬ 
portional to (T-p and the local x density pj^: 


Xoiafpp^. 


(24) 


In the case where the CDM consists only of neutralinos. 


Px = PCDM, 


(25) 


where pcDM is the local CDM density. So, the LUX 
experiment bound on A can be directly translated into a 
bound on provided that pcDM can be estimated - see 
e.g. Ref. |l03j| . However, when the x’s constitute only a 
fraction of the total CDM, we can write Eq. (IMl) as 

A oc cr|p PCDM- (26) 

PCDM 

Following the authors of Ref. |in4l | , we can then use the 
scaling ansatz 


Px/PCDM = AI^/Dcdm, 


(27) 


which gives 


A oc (t|p ^ PCDM. (28) 

“CDM 

So, the LUX experiment bound on A is now translated 
into a bound on the ‘rescaled’ SI neutralino-proton elastic 
cross section Ccr|p, where ^ = Hj^/HcDM- 


III. THE HYPERBOLIC BRANCH/FOCUS 
POINT REGION 

A detailed discussion of the characteristics of the 
HB/FP region of the CMSSM is given in Refs. [H, 

The classification of the various solutions of the radiative 
EWSB condition is based on the expansion of in terms 
of the soft SUSY breaking parameters of the CMSSM in¬ 
cluded in Eq. O- Indeed, using fitting techniques, we 
can verify the following formula 

-\- M '^12 ~ cottiq -\- €112X1^^2 + caAq -|- camXqMiJ 2 ^ 

(29) 

where the coefficients co,Ci/ 2 ,ca, and cam depend ba¬ 
sically on tan/3 and the masses of the fermions of the 
third generation. These coefhcients are computed at the 
scale Msusy in Eq. ([3]) and, therefore, inherit a mild de¬ 
pendence on the SUSY spectrum too. From Eq. (1331) . 
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TABLE I: The c’s in Eq. (1291) for the four cases of Table ITIl 


tan /3 

Co 

Sco 

Cl/2 

CA 

Cam 

40 

-0.0921 

0.014 

0.789 

0.107 

-0.269 

45 

-0.0775 

0.0148 

0.825 

0.1016 

-0.260 

48 

-0.0686 

0.0151 

0.845 

0.0981 

-0.253 

50 

-0.0624 

0.0265 

0.859 

0.0953 

-0.247 


we can easily infer that the SUSY breaking parameters 
are bounded above for fixed /r, when the quadratic form 
in the right-hand side of this equation is positive defi¬ 
nite. This is the so-called ellipsoidal branch (EB). On 
the contrary, in the hyperbolic branch (HB) region 
Co is negative and, consequently, mo can become very 
large together with a combination of Aq and Mi ^2 with 
all the other parameters being fixed. Near the boundary 
between the EB and HB regions, the coefficient cq is very 
close to zero and, thus, mo can become very large with 
all the other parameters fixed. Moreover, there is a re¬ 
gion where the soft SUSY breaking mass-squared of 
H 2 becomes independent of the asymptotic value of the 
parameter rriQ. This is called the focus point (FP) region 

In the HB region, the radiative EWSB admits three 
types of solutions [l^, : 

(i) Focal Points (HB/FP): They lie near the boundary 

between the EB and the HB regions, where co — 0 
and, thus, p? is practically independent of m^. As 
a consequence, m^ can become very large, while 
pf and all the other parameters remain fixed. We 
should note that, in the large tan/3 regime, we ob¬ 
tain p^ + M‘^12 ~ focal points, 

where p'^ is essentially independent of mg, coincide 
with the focus points, where is almost mg- 

independent. Since, in the present model, tan/3 
is large, we will not distinguish focal from focus 
points and we will use the same abbreviation (FP) 
for both. 

(ii) Focal Curves (HB/FC): Along these curves, two of 
the soft parameters can acquire large values, while 
p^ and all the other parameters remain fixed. 

(hi) Focal Surfaces (HB/FS): Here we can have three 
soft parameters with large values, while p'^ and the 
other parameters remain fixed. 

To get an idea of how our solutions in Sec.|lV]are clas¬ 
sified into these categories, we display, in Table HI the val¬ 
ues of the coefficients c in Eq. (IMl) for the four representa¬ 
tive cases of Table HTl corresponding to tan (3 = 40,45,48, 
and 50 - see Sec. IIVDI We also list the relevant 1 — cr 
uncertainly Scq of cq, which is practically the deviation 
between the lowest cq obtained by varying mt{mt) in the 
1 — tj range in Eq. 0 and the value of cq correspond¬ 
ing to the central value of mt(mt) - possible variation of 
mb(Mz) in its 1 — cr range in Eq. ® does not change cq 


significantly. As it turns out, the minimal cq corresponds 
to the upper limit of mt(mt) in Eq. ([7]). Taking then the 
ratio |co|/i5co, we see that the displacement of cq from 
zero ranges from 2.4 — a (for tan ,5 = 50) to 6.6 — cr (for 
tan/3 = 40). In particular, for tan/3 > 45, this displace¬ 
ment does not exceed 5 — cr and, thus, a real deviation of 
the coefficient cq from zero cannot be established. As a 
consequence, we can say that the corresponding param¬ 
eters lie within the HB/FP region. On the contrary, for 
tan/3 < 45, the parameters belong to the HB/FC and 
HB/FS region. 


IV. RESTRICTIONS ON THE SUSY 
PARAMETERS 

Imposing the requirements described in Sec. |TT] ~ with 
the exception of the one of Sec. Ill A 5I -. we can restrict 
the parameters of our model. Following our approach in 
Refs. [ 4 ^, we use as free parameters of the model 
the ones in Eq. O- The Yukawa coupling constant ra¬ 
tios ht/hr and hb/hr are then fixed by using the data 
of Sec. [irXTl These ratios must satisfy the YQU con¬ 
ditions in Eq. @ for natural values of the parameters 
ai,a 2 , and p - see Sec. El We restrict ourselves to the 
range 40 < tan/3 < 50 since, below tan/3 = 40, the ratios 
hm/hn {m,n = t,b,T) tend to require unnatural values 
of ai,a 2 , and p and, above tan/3 = 50, the numerical 
calculations for the soft SUSY masses become quite un¬ 
stable. As a characteristic value in our allowed parameter 
space, we take tan/3 = 48, which balances well enough 
between maintaining natural values for the hm/hn’s and 
satisfying the various requirements of Sec. [TTl 

We concentrate on the p > 0 case, given that ^ < 0 
worsens the violation of Eq. CZiD, and scan the region 
—30 < AofMi /2 < 30. We will hereafter use the central 
values of the SM parameters Mj, mb{Mz), mr{Mz)^ and 
as{Mz) given in Sec. Ill All We find that the only con¬ 
straints which play a role are the lower bound on mb in 
Eq. (fTUl) . the bounds on m^± and nig in Eqs. (fTTl) and 
m, the CDM bound in Eq. m, and the bound from 
the LUX experiment - see Sec. Ill B 21 In the param¬ 
eter space allowed by these requirements, all the other 
restrictions of Sec. Ill Al are automatically satisfied with 
the exception of the lower bound on Sa^ in Eq. (|17bl) . 
This bound is not imposed here as a strict constraint on 
the parameters of the model for the reasons explained 
in Sec. Ill A 51 We only discuss the level at which the 
requirement in Eq. (I17al) is satisfied in the parametric 
region allowed by all the other constraints. 

We first present, in Sec. IIV Al the restrictions in the 
Mi /2 ~ plane for low values of Ag/Mj ^/2 then, in 
Sec. IIVBI we reveal the role of the bound in Eq. (Ha) by 
considering the tan /3 = 48 case and ignoring the restric¬ 
tions of Sec. Ill B 21 In Secs. IIV Cl and IIVDI we present 
regions for various tan/3’s consistent with all the require¬ 
ments of Sec. m but the one of Sec. Ill A 51 
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tan|3=48 



FIG. 2: Allowed (shaded) regions in the M 1/2 — rno plane for 
tan/3 = 48 and different indicated in the plot. The 

red lines correspond to rrih in Eq. m- 


A. Restrictions in the M1/2 — mo Plane 

The interplay of the various requirements of Sec.HIlcan 
be easily understood from Fig. [2l where we present the 
(shaded) strips in the M 1 /2 — W q plane allowed by the 
restrictions of Secs. nTXTirrr^irT^nTXdi and irfRii 
for tan/3 = 48 and several values of A 0 /M 1/2 indicated 
in the graph. The upper boundary along each of these 
allowed strips arises from the limit on in Eq. (HU. 
Note that this limit is more restrictive than the limit for 
triggering radiative EWSB - cf. Ref. |l05|| . The lower 
boundary along each strip is given by the upper limit on 
ili^sph^ in Eq. m- On the other hand, the lower limit on 
m/i in Eq. (jH]) causes the termination of the strips at low 
values of m-o and M 1 / 2 , whereas their termination at high 
values of mg is put in by hand in order to avoid shifting 
the SUSY masses to very large values. The red lines 
indicate solutions with nih = 125.5 GeV - see Eq. (nnj. 

From this figure, we easily see the main features of the 
HB/FP region: mg spans a huge range (4 — 15) TeV, 
whereas M 1/2 [m] remains relatively low (1 — 6 ) TeV 
[(0.1 — 1) TeV] thanks to the structure of the radiative 
EWSB in this region analyzed in Sec. lIIIl We observe also 
that as Ao/Mi /2 increases from —2 to 2 the allowed strip 
moves to larger Mi/ 2 ’s and becomes less steep, while, for 
AqIMi ^2 > 2 , we have exactly the opposite behavior: as 
AqIMiI 2 increases, the allowed strip moves to smaller 
values of M 1/2 and, at the same time, becomes steeper. 


B. Restrictions from the Bound on mg 

As IA 0 /M 1 / 2 I increases above 10, the restriction on 
mg in Eq. m comes into play. Its importance can be 
easily inferred from Fig. [31 where we plot mg against 
Ao/Mi /2 for mh as in Eq. (fTUl) and tan/3 = 48 - let us 


tanp=48 



FIG. 3: Restrictions in the A 0 /M 1/2 — mg plane from Eq. (Hill 
(solid line) and Eq. (I19II (dashed line) for mg in Eq. (IIUII and 
tan /3 = 48. The region excluded by Eq. m is cross hatched. 


note that varying tan /3 does not significantly affect this 
figure. The solid [dashed] line corresponds to the limit 
on mj^± [flnsp^^] in Eq- (flTI) [Eq. (IT^ ]. The dotted line 
denotes the lower limit on mg in Eq. (1121) excluding the 
cross hatched region. Therefore, the allowed region (still 
without considering the LUX data - see Sec. Ill B 2|) is 
the region bounded by the three aforementioned lines. 
We see clearly that the bound on mg cuts off the largest 
values of lAg/Mi/ 2 |. Specifically, near the lower bound 
on m^± given in Eq. (fTTl) . Ag/Mi /2 is limited in the 
range —10.8 ;< Ag/Mi /2 14.4, while, near the up¬ 
per bound on in Eq. (HU, it is limited in the 

range —11.7 < Ag/Mi /2 < 15.5. In order to understand 
this behavior, let us recall that mg depends on the ratio 
/m?, where Xt = At — fj, cot /3 with At being the soft 
trilinear scalar coupling constant for the t-squark. For 
mg and tan /3 fixed, we see that larger absolute values 
of Ag/Mi /2 require larger m^’s too and, since mi de¬ 
pends largely on mg, the latter must grow larger as well 
together with the left-handed squark soft SUSY break¬ 
ing masses Mq^ . These soft masses increase as the soft 
gluino mass M 3 decreases - see the relevant renormal¬ 
ization group equation in Ref. Q - and, thus, it is quite 
obvious that, in order to achieve larger values of Mq^, 
we need smaller Mg’s, which also suggests smaller values 
for the common asymptotic gaugino mass Alij 2 - There¬ 
fore, small Mi/ 2 's and large J^cJ’s, which lead to big 
jAg/Mi/ 2 |’s, yield small mg’s, which is exactly what one 
deduces from Fig. [3] 


C. Restrictions from the LUX Experiment 

Considering the cross section tr-p, we can not only 
probe the detectability of the LSP, but also obtain fur¬ 
ther restrictions on the parameters of our model. This 




















is because of the rather enhanced cr|p’s obtained in the 
HB/FP region, as explained in Sec. Ill B 21 In the com¬ 
putation of cr~p, we adopt the central values of the /^g’s 
in Eqs. (I22al) - (l22cl) and fix ruh to its value in Eq. 

The results are presented in Fig. lU where we depict the 
allowed (bounded) regions in the tolsp — plane for 
tan /? = 40 (upper panels) tan /3 = 48 (middle panels) and 
tan/3 = 50 (lower panels). We also draw with gray solid 
and dashed lines the projected sensi tiviti es of the 
XENONIT [Toi] and LUX-ZEPLIN [To^ experiments, 
respectively. The panels in the left [right] column corre¬ 
spond to ^o/Mi /2 < 0 [Ao/Mi /2 > 0]. We see that, for 
each tan /3, the allowed regions in the left and right panel 
almost coincide and, thus, we display them separately to 
avoid confusion. The numbers on the various points of 
each boundary line indicate the corresponding values of 
AojMii2. 

The allowed regions are bounded by five different types 
of black lines for which we adopt the following conven¬ 
tions: 

• on the solid line the bound on m^± coming from 
Eq. (fTTl) is saturated; 

• on the double dot-dashed line the limit on rrig in 
Eq. (fT^ is saturated; 

• on the dashed line the bound on from 

Eq. (fTOl) is saturated; 

• the dotted line depicts the bound on arising 
from the LUX data taken from Ref. [^ : 

• the dash-dotted line represents the lowest possible 

in each case. 

As can be understood from the description of the various 
lines above, we have ^ = 1 along the dashed line. As 
regards the minimal this reaches 0.013 independently 
of tan/3. 

Notice that there are no significant differences between 
the various tan /3’s as regards the allowed ranges of mnsp 
and AqIMi/ 2 - From all these graphs, we deduce that, 
as we move towards the Ulsp^^ = 0.125 line, tolsp and 
mj^± grow larger, while the allowed range of A 0 /M 1/2 
becomes smaller. The largest IA 0 /M 1 / 2 I is located at the 
junction point of the double dot-dashed and dotted lines 
- with the exception of the allowed region for tan /3 = 50 
and Ao/Mi /2 > 0, where the largest IA 0 /M 1 / 2 I is found 
at the intersection of the solid and double dot-dashed 
lines. The smallest mnsp can be found in the upper left 
corner of the allowed regions, whereas the largest tolsp is 
found on the dashed line. The allowed ranges of A 0 /M 1/2 
and Tolsp can be summarized as follows: 

• For tan/3 = 40, we find —12.4 < A 0 /M 1/2 ^ 16.28 
and 92 < mnsp/GeV < 1084.8. 

• For tan/3 = 48, we find —12.8 < A 0 /M 1/2 ^ 15.8 
and 92 < mnsp/GeV < 1084.2. 


• For tan/3 = 50, we find —13 < A 0 /M 1/2 ^ 15.35 
and 91.9 < mnsp/GeV < 1088. 

We noticed that, as A 0 /M 1/2 approaches 2 from above, 
the value of the product |A^i,ip|A^i, 3 p decreases. How¬ 
ever, as Aq/Mi /2 decreases below 2, |Afi,ip|A^i, 3 p grows. 
The growth of this product persists even when A^jM^j^ 
becomes negative. As a consequence - see Eq. (1^ ~ 'fcr|p 
acquires its minimal value at A 0 /M 1/2 = 2. The over¬ 
all minimum of ^cr-p, for each value of tan/3, is acquired 
at the lowest left corner of the corresponding allowed re¬ 
gion in the right column of Fig. H] The smallest of these 
overall minima, for tan/3 > 40, is 

^crSI > ^ j^q- 12 4g ^ ^q-12 p^^^ ^ 

and it is obtained at tan /3 = 40 and ttilsp — 101.7 GeV. 
The bound in parenthesis is derived by allowing the /^n’s 
to vary within their 1 — cr intervals in Eqs. (I22a|) . (I22bl) . 
and H22cD . As shown in the plots, the obtained values 
of within the reach of forthcoming experiments 

like XENONIT and LUX-ZEPLIN with planned sensitiv¬ 
ity from 10“"^® to 10“"^^ cm^ for the ttilsp range under 
consideration - recall that 1 pb = 10“^® cm^. 

D. The Overall Allowed Parameter Space 

In order to have a more spherical view of the allowed 
parameter space of our model, we present in Fig. [5| the 
allowed regions in the ttia — A 0 /M 1/2 plane for tan/3 = 
40,48, and 50 enclosed by blue, red, and green lines, re¬ 
spectively. As one can see from Fig. |4| and the allowed 
ranges of thlsp presented in Sec. IIV Cl the mass of the 
LSP (and M 1 / 2 ) is confined in approximately the same 
range independently of tan/3. Therefore, we opt to use 
mA as variable in the horizontal axis, so that the al¬ 
lowed regions corresponding to different tan/3’s are well 
distinguishable. We observe that, although the range 
of Tolsp is practically tan/3-independent, mA increases 
as tan,5 decreases. For each value of tan/3, there are 
three different boundary lines corresponding to different 
restrictions. On the solid and dashed line, the bounds on 
in Eq. dD and on Ulsp^^ in Eq. m are saturated, 
whereas the restriction from the LUX data on yields 
the dotted boundary line. It is impressive that the LUX 
data provide such a strong constraint on the model pa¬ 
rameters, which overshadows all other constraints for ap¬ 
proximately Ao/A/ 1/2 ^ — 3 and Ao/A'/i /2 ^ 5. Finally, 
the boundary lines from the limit on mg in Eq. (1121) . al¬ 
though hardly visible in this plot, provide the maximal 
and minimal Ao/Mi/ 2 ’s. Note that the allowed regions 
are obviously symmetric about A 0 /M 1/2 — 2.5. Also, 
we find that fjL remains almost constant ~ 100 ± 20 GeV 
on the solid lines from Eq. (HU, while it reaches about 
1 TeV when the bound in Eq. HU) is saturated. As re¬ 
gards <3a^, the acquired values are well below the lower 
limit in Eq. (I17bl) . Specifically, in the allowed regions of 
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tanp=40 tanp=40 




tanp=48 



tanp=48 



tanp=50 



itilsp (GeV) 


tanp=50 



itilsp (GeV) 


FIG. 4: Allowed (bounded) regions in the ttilsp — plane for rrih given by Eq. (IIOII and tan {3 = 40, 48, and 50. The left 
[right] panels correspond to < 0 [A 0 /M 1/2 > 0] and the values of AqIMi /2 at the various points of the boundary lines 

are indicated. The black solid, double dot-dashed, and dashed lines correspond to the bounds from Eqs. (dll), (HU), and (O, 
respectively. The black dotted lines arise from the LUX data, whereas the black dot-dashed lines give the lowest possible 
in each case. The planned sensitivity limits of XENONIT and LUX-ZEPLIN are also depicted by a gray solid and dashed line, 
respectively. 
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FIG. 5: Allowed regions in the rriA — A 0 IM 1/2 plane with m/, 
given by Eq. (nsj and for tan /3 = 40,48, and 50, bounded by 
blue, red, and green lines, respectively. The solid, dashed, and 
dotted boundary lines come from the constraints in Eqs. (HU, 
(HU, and the LUX data, respectively. The restriction in 
Eq. m limits the range of A 0 IM 1 / 2 , but the corresponding 
lines are hardly visible. 

Fig.[^ we obtain <5o^ ~ (0.04 — 0.27) x 10“^°. Therefore, 
Eq. (jl7al) is satisfied only at the level of 2.83 to 2.86 — a. 

The values of the input and some output parame¬ 
ters, the mass spectra and some low energy observ¬ 
ables of the model are listed in Table El for four char¬ 
acteristic points of the allowed parameter space with 
rrih = 125.5 GeV. The various masses of the SUSY parti¬ 
cles (gauginos/higgsinosx,^X2, X3,X4,xf,X2> 9, squarks 
ii, ^2, &i, &2, ul, Ur, dL, dn, and sleptons n, f2, f'r, sl, 
en) and the Higgs particles {h, H, A) are given in 
TeV - note that we consider the hrst two generations of 
squarks and sleptons as degenerate. We chose the values 
of the input parameters so as to assure that is 

not far from its central value in Eq. (m. The relatively 
low A^i’s and Ay^’s obtained and the sizable higgsino 
fraction |iVi^3p + assist us to achieve this. Note 

that, in all the cases, the lightest neutralino mass turns 
out to be close to the corresponding /r. We also include 
an estimate for the EWSB fine-tuning parameter Aew “ 
see Sec. I VII 

We observe that the low energy B-physics observables 
satisfy the relevant constraints in Sec. Ill whereas the 

(5a^’s are far below the lower limit in Eq. (I17bl) . Also, 
the extracte d ca n be probed by the forthcoming 

experiments |l06L llOTf . For the spin-dependent x — p 
scattering cross section , we adopt the central values 
for the hadronic inputs in Ref. [5^. As can be easily 
deduced from the displayed values, in our model lies 
well below the sensitivity of IceCube |l08j| (assuming y — 
X annihilation into W^W~) and the expected limit from 
the large DMTPC detector [9^. Therefore, the LSPs 


predicted by our model can be detectable only by the 
future experiments which will release data on 

Comparing the results of Table HIl with the correspond¬ 
ing ones in the x — U coannihilation region - see Table II 
of Refs, 113 we may appreciate the different fea¬ 
tures of the HB/FP solutions presented here. First of 
all, mo acquires considerably larger values here, while p 
remains quite small. The Higgs bosons H, and A 
acquire larger masses and the whole sparticle spectrum, 
with the exception of the neutralinos and charginos, be¬ 
comes heavier. The gluino mass is heavier for lower 
tan /3’s, in contrast with what happens in the x — ti coan¬ 
nihilation region, where we observe the opposite behav¬ 
ior. Here, the ratio hi,/hr is smaller, ht/hi, is closer to 
2, while ht/hr is even closer to unity. The restrictions 
on the low energy observables are all well satisfied, ex¬ 
cept the one on Sa^, which becomes even smaller than 
in the x — ti coannihilation region. Note that the lat¬ 
ter region is tightly constrained by BR(Bs —^ pAp~), 
which is well suppressed in the HB/FP region. As re¬ 
gards this becomes compatible with Eq. (ITO)) in 

the present case thanks to x/X2 ~ Xi" coannihilations, 
which are activated because of the low A^±’s and A^^ ’s. 
Note that, in the x ~ U coannihilation region, the rele¬ 
vant mass splittings are Af , = (ruf^ — tolsp)/™lsp and 
Ah = {mn — 2TOLSp)/2m.LSP- Finally, here we obtain a 
large higgsino fraction of the LSP, which confines mnsp 
close to pL and as a bonus ensu r es al ^’s accessible to the 
forthcoming experiments [IO6L Il07j |. On the contrary, 
within the x — ti coannihilation region, x is an almost 
pure bino, mnsp — M 112 / 2 , and cr|p is well below the 
sensitivity of any planned experiment. 

V. DEVIATION FROM YUKAWA 
UNIFICATION 

The ranges of the ratios of the asymptotic third genera¬ 
tion Yukawa coupling constants in the allowed parameter 
space of the model in the range 40 < tan/3 < 50, which 
we consider here, are the following: 

1 ^ ^ ^ 1-5, 0.75 <^< 0.79, and 1-2 < ^ < 2. 

ti-T fij- hi) 

(31) 

It is easy to see that, although YU is violated, these ra¬ 
tios turn out to be quite close to unity. However, for 
tan/3 < 40, the ratios ht/hr and ht/hb become large and 
the term ‘Yukawa quasi-unification conditions’ cannot be 
justified for the constraints given in Eq. ©■ So, we re¬ 
strict ourselves to tan /3 > 40. 

We will now give a specific example of how such values 
for the ratios of the Yukawa coupling constants as the 
ones in Eq. m can be achieved in a natural way. We 
consider the third case in Table HIl with tan/3 = 48, where 
ht/hr = 1.107 and h^/hr = 0.763, and solve Eq. @ with 
respect to the complex parameters ai, 02, and the real 
and positive parameter p. Specifically, we first find pairs 
of values for p and ai which satisfy the YQU condition 
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TABLE II: Input/output parameters, sparticle and Higgs masses, and low energy observables in four cases in the HB/FP. 


Input parameters 

tan /3 

40 

45 

48 

50 

Ao/Mi/2 

2 

0 

-1.5 

-3 

Mi/ 2 /TeV 

3.763 

2.945 

2.161 

1612.4 

mo/TeV 

9.603 

8.821 

9.231 

9.561 

Output parameters 

ht/hr{MGlJT) 

1.474 

1.237 

1.107 

1.027 

hb/hr{MGVT) 

0.756 

0.758 

0.763 

0.774 

ht/hb{MGVT) 

1.949 

1.631 

1.45 

1.326 

p/TeV 

1.018 

1.01 

0.928 

0.736 


0.002 

0.004 

0.014 

0.038 

Ax. 

0.004 

0.006 

0.018 

0.045 

Aew 

244.6 

237.9 

216.2 

130.9 

|Ni,3|^ + \Nl,4f 

0.997 

0.992 

0.836 

0.54 

Sparticle and Higgs boson masses in TeV 

X,X2 

1.043,1.047 

1.026,1.032 

0.935,0.952 

0.723,0.756 

X3,X4 

1.792,3.381 

1.397,2.651 

1.034,1.959 

0.788,1.473 


1.046,3.381 

1.030,2.651 

0.949,1.959 

0.751,1.473 

g 

8.147 

6.512 

4.951 

3.799 

tl, t2 

8.113,9.516 

6.900, 7.984 

6.309, 7.270 

6.049, 6.900 

bl, 1)2 

9.515,10.23 

7.987, 8.560 

7.267, 7.887 

6.897, 7.512 

UL^ UR 

11.777,11.553 

10.357, 10.197 

10.100,10.004 

10.082,10.022 

dL, dit 

11.777,11.524 

10.357, 10.176 

10.100,9.992 

10.082,10.014 

n, f 2 

7.927,9.104 

6.921,8.136 

6.749, 8.202 

6.620, 8.292 

Vr 

9.103 

8.135 

8.201 

8.291 

Ci, e.R 

9.936,9.718 

9.050,8.900 

9.359,9.276 

9.640,9.589 


9.935 

9.049 

9.359 

9.639 

h,H 

0.1255,6.56 

0.1255,4.820 

0.1255,3.67 

0.1255,2.369 

H^,A 

6.56,6.56 

4.820,4.820 

3.671,3.67 

2.370,2.371 

Low energy observables 

lO^BR {b 57 ) 

3.31 

3.30 

3.30 

3.33 

10®BR (i?s ^ p+p“) 

3.04 

3.03 

3.02 

2.97 

R(i3„ Tu) 

0.998 

0.995 

0.991 

0.976 

W^°6a„ 

0.135 

0.2 

0.227 

0.237 


0.11 

0.11 

0.11 

0.11 

CT|p/10“®pb 

0.28 

0.81 

7.75 

13.29 

ct|p/ 10~^pb 

2.55 

7.87 

77.08 

227.8 


for hb/hr and, then, for each one of these pairs, we find 
0 : 2 ’s which satisfy the equation for ht/hr- Note that, 
while hb/hr depends only on the value of the product 
pai, the ratio ht/hr depends on each of the parameters 
ai, a 2 , and p separately. This is the reason why the 
solutions in the ai complex plane lie on a set of similar 
curves corresponding to different values of p, whereas the 
values of 0:2 do not follow any specific pattern in the 
corresponding complex plane - see below. We found that 
solutions exist only for low values of p (up to about 0 . 6 ) 
and also for large values of this parameter (p > 2.4). 

We present several of these solutions in the ai and 0:2 
complex planes in Fig. [51 but only for the lower values of 
p, which are considered more natural. In the left panel 


of this figure, we see that, for any given value of p, the 
data clearly lie on a specific curve and that the curves 
for different values of p are similar to each other. On 
the other hand, in the right panel, the data are more 
complicatedly distributed. Note that, every pair of values 
of ai and p corresponds to more than one value of a 2 - 
From Fig. we can deduce that, for the specific example 
considered and for small values of p, the ranges of the 
values of all the parameters are the following: 

0-3 < p <0.5, (32a) 

-0.45 < Re(ai) < 0.31, 0.95 < Im(ai) < 2.19, (32b) 
-1.39 < Re(a 2 ) < 0.83, 0.01 < Im(a 2 ) < 1.21. (32c) 
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FIG. 6: The complex parameters oi and a 2 for various real and positive values of p indicated on the graphs for the case in the third 
column of Table HD with tan /3 = 48. 


Note also that, for the larger values of p, we get smaller 
values for ai and a^- So, we see that the ratios of the 
Yukawa coupling constants in this example can be easily 
obtained by natural choices of p, ai, and 02 - 

We hnd that, for all the possible values of the ratios of 
the third generation Yukawa coupling constants encoun¬ 
tered in our investigation, the picture is quite similar. So, 
we conclude that these ratios can be readily obtained by 
a multitude of natural choices of the parameters p, oi, 
and 0.1 everywhere in the allowed parameter space of the 
model which we considered here. 


VI. NATURALNESS OF THE EWSB 


One of the most important reasons for introducing 
SUSY was that it could provide a solution to the hierar¬ 
chy problem. However, the fact that, in our model, toq 
turns out to lie in the multi-TeV range generates doubts 
about the naturalness of the radiative EWSB, since it 
leads to the reappearance of a small fine-tuning problem. 
This is the so-called little hierarchy problem, which is still 
a much debated issue. 

To quantify somehow the naturalness of our model 
with respect to this issue, we focus on the EWSB con¬ 
dition relating Mz to p and tan ,5. This condition is 
obtained by minimizing the tree-level renormalization- 
group improved scalar potential for Hi and H 2 and reads 
as follows 0 


1 A ^2 /3 2 

tan2/3-l 


(33) 


where ruHi is the soft SUSY breaking mass of Hi and 
we neglect possible loop corrections, which are anyway 


minimized thanks to the choice of the optimal scale in 
Eq. ([3]). In the case where the value of the fraction in 
the right-hand side of Eq. (ESI is quite large, a large 
value of is also needed. However, this requires large 
values of mo and M 1/2 and heavy sparticle spectrum and, 
therefore, introduces a certain amount of fine-tuning. To 
measure this tuning, we introduce the EWSB fine-tuning 
parameter 

AEw^max(^), (34) 

where the Q’s {i = p,Hi,H 2 ) are - see Eq. (l33l) - : 


[C^,Ch,.Ch2) 


m 


Hi 


tan^ 13—1 


mf^^ tan^ (3 \ 

tan^ /3 — 1 J 

(35) 


In most of the parameter space explored, Aew is domi¬ 
nated by the term 

Focusing on the values of the parameters which ensure 
= 0.125, we present, in the left panel of Fig. [71 
Mi /2 (solid line) and Aew (dashed line) as functions of 
wlsp for tan (3 = 48, mu = 125.5 GeV, and negative val¬ 
ues of AoIMi /2 indicated in the graph - cf. the left panel 
in Fig.|3]for tan/3 = 48. We clearly see that the required 
EWSB fine-tuning parameter Aew is almost constant 
and ~ 200. This result is also valid for the other two 
values of tan/3 in Fig. U) Indeed, for the largest value of 
IA 0 /M 1 / 2 I on the dashed line in each of the six panels of 
Fig. m the derived Aew is as follows: 


• For tan/3 = 40 and A 0 /M 1/2 = —1.6 [A 0 /M 1/2 = 
5.7], we find Aew = 191.311 [Aew = 187.379]. 

• For tan/3 = 48 and A 0 /M 1/2 = —1.9 [A 0 /M 1/2 = 
5.6], we find Aew = 188.064 [Aew = 187.724]. 
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FIG. 7 : Mi^2 Aew functions of rriLSP for — 0 . 125 , TUh given by Eq. {Toll and various Ao/Mi/2’s indicated on the curves 

for the HB/FP (left panel) and the x ~ "ri coannihilation (right panel) region of the model. 


• For tan/3 = 50 and = —1-8 [A0/M1/2 = 

5.2], we find Aew = 204.702 [Aew = 201.228]. 

The value of Aew can be even smaller for lower ttilsp’s 
and riLSP^^’s. However, we believe that the values pre¬ 
sented above are more interesting, especially if we wish 
to compare the Aew’s found in the HB/FP region with 
the ones in the x — ti coannihilation region analyzed in 
Ref. [5^. To this end, we plot also, in the right panel 
of Fig. 0 Mi /2 (solid line) and Aew (dashed line) as 
functions of tolsp in the x ~ coannihilation region of 
the same model for tan/3 = 48, nu = 125.5 GeV, and 
AqIMii 2 < 0 - see Fig. 3 of Ref. [s^. It is evident that, 
in this coannihilation region, the required values of ttilsp 
are somewhat larger and the values of Aew can become 
about a factor of ten larger (Aew > 1000). We can, thus, 
easily appreciate the amelioration regarding the EWSB 
fine-tuning that we achieve working in the HB/FP region. 

VII. CONCLUSIONS 

We investigated the compatibility of the generalized 
asymptotic Yukawa coupling constant quasi-unification, 
which yields acceptable masses for the fermions of the 
third family, with the HB/FP region of the CMSSM for 
/r > 0 and 40 < tan /3 < 50. We imposed phenomenolog¬ 
ical constraints originating from the mass of the lightest 
neutral CP-even Higgs boson, the lower bounds on the 
masses of the sparticles, and H-physics. We also con¬ 
sidered cosmological constraints coming from the CDM 
relic abundance in the universe and the LUX data on 
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